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Q ! ABSTRACT 

^ ■ Mapping CMB polarization is an essential ingredient of current cosmological re- 

search. Particularly challenging is the measurement of an extremely weak B-mode 
polarization that can potentially yield unique insight on inflation. Achieving this ob- 
jective requires very precise measurements of the secondary polarization components 
on both large and small angular scales. Scattering of the CMB in galaxy clusters in- 
I duces several polarization effects whose measurements can probe cluster properties. 

. Perhaps more important are levels of the statistical polarization signals from the popu- 

O^l ' lation of clusters. Power spectra of five of these polarization components are calculated 

^SJ . and compared with the primary polarization spectra. These spectra peak at multipoles 

' £ ^ 3000, and attain levels that are unlikely to appreciably contaminate the primordial 

0^ . polarization signals. 
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1 INTRODUCTION 

Primordial seed perturbations in the densities of radiation and matter, and the excitation of gravitational waves during an 
epoch of cosmological inflation, left their imprints on the cosmic microwave background (CMB). Generic models of inflation 
predict gaussian density fluctuations, and therefore all the statistical information characterizing the primordial CMB is 
expected to be encoded in the two-point correlation function defined on the celestial sphere, or equivalently in its harmonic 
transform - the angular power spectrum. Various processes that took place at later times when the evolution on small scales 
was nonlinear, can induce non-gaussian anisotropy. An important example is the anisotropy induced by Compton scattering 
of the CMB by hot gas in galaxy clusters - the SZ effect. Indeed, for this effect higher order spectra, e.g. the tri-spectra, 
may contain additional valuable information. Nonetheless, a wealth of information can be gleaned from the angular power 
spectrum itself. This work focuses on power spectra of CMB polarization induced by scattering in clusters. 

While the primary temperature anisotropy, velocity-gradient-induced E-mode polarization, and gravitational-wave- induced 
B-mode polarization dominate the lowest and intermediate multipoles of the power spectra, the latter is damped on scales 
smaller than the horizon at recombination, and the former are exponentially suppressed beyond I « 1000 due to photon 
diffusion damping (Silk 1968). The latter process took place during recombination, when the radiation decoupled from the 
baryons for the first time. On these few arcminute scales (determined by the thickness of the last scattering shell) , secondary 
signals induced by clusters of galaxies peak at multipoles, I, characteristic of galaxy clusters (/ « 2000 — 3000). This secondary 
anisotropy gauges the formation and evolution of clusters; as such it depends on different combinations of the cosmological 
parameters than the primary CMB anisotropy. Its steep dependence on the basic quantities that determine the evolution of 
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the large scale structure (LSS) - such as the critical density for collapse and mass variance - makes it particularly valuable for 
removing some of the inherent degeneracies in the cosmological parameters. 

The generation of temperature anisotropy by clusters can be either through Icnsing by their gravitational field, which is 
clearly dominated by the dark sector of the cluster total energy (i.e. dark matter, and possibly also - indirectly - dark energy), 
or by Compton scattering of the radiation by intracluster (IC) gas, the main compoent of the baryonic sector. Cluster-induced 
Icnsing of the temperature and polarization anisotropy was the subject of several recent studies. It merely redistributes 
the anisotropy on the sky while conserving the total power. Lensing couples fluctuations of temperature or polarization on 
cluster scales to the lensing deflection angle, which can be directly related to the transverse gradient of the cluster projected 
gravitational potential. Therefore, lensing smoothes features in the power spectrum on scales comparable to the characteristic 
deflection angle. Lensing of the CMB by clusters also converts the polarization modes from E to B, which is indeed part of 
the full lensing signature of the LSS. 

In this paper we focus on the statistical polarization signals induced collectively by scattering in clusters. Temperature 
anisotropy is induced when the CMB is scattered by moving electrons; fully random motions give rise to the thermal SZ 
effect, and additionally also to the kinematic SZ effect when the cluster has a flnite radial velocity in the CMB frame. 
The kinematic effect is typically an order of magnitude smaller than the thermal component. Additional information can 
in principle be extracted from the polarization state of the radiation, which probes other combinations of cosmological and 
cluster parameters. 

Various polarization effects were studied in the original work of Sunyaev & Zeldovich (1980), and further elaborated 
upon by Sazonov & Sunyaev (1999) and e.g. Audit & Simmons (1999). Conversion of the primordial E-mode polarization to 
B-mode by the LSS has been considered by e.g. Seljak & Zaldarriaga (2000), Lewis & King (2006) and by Hu, DeDeo & Vale 
(2007). Gravitational lensing of the primary CMB mixes the two E & B polarization which makes their separation from the 
primordial polarization extremely challenging, necessitating the use of higher order statistics to de-lens the sky. While these 
delensing methods work well for the LSS (because they were optimized for this case), similar methods should be applied to 
remove cluster signals; however, this need may be circumvented if the scattering-induced polarization signals are very weak, 
as seems to be the case based on results presented in this paper. 

Cluster-produced polarization signals are typically much smaller than the temperatore anisotropy by orders of magnitude. 
The only two such polarization components discussed so far in the context of their statistical imprint on the polarization power 
spectrum are those resulting from the primordial quadrupole and the quadrupole anisotropy associated with the transversal 
(Doppler) component of the bulk motion of the cluster (Cooray, Baumann & Sigurdson 2005). The latter is second order 
in the cluster transverse velocity, and linear in the optical depth. Similarly, we expect the quadrupole induced by tensor 
perturbations to generate polarization, but this effect is much smaller, as implied from the current upper limit, 0.22 (95% 
CL), on the tensor-to-scalar ratio (Hinshaw et al. 2008) 

Two other relevant polarization components whose statistical properties were not studied arise from double scattering 
in the same cluster (Sunyaev & Zeldovich 1980, Sazonov & Sunyaev 1999). First scattering induces temperature anisotropy 
either by the thermal or kinematic SZ effect. If this temperature anisotropy contains a local quadrupole moment, polarization 
is induced upon second scattering. For a typical Thomson optical depth of IC gas, r ~ 0.01, these r^-depenedent components 
are clearly very weak, with the thermal effect sourcing the largest of the two. We ignore here the minute signals due to 
radiation scattered during aspherical collapse of protoclusters, scattering in a moving cluster in which the radiation develops 
anisotropy due to Icnsing (Gibilisco 1997), or polarization produced by IC magnetic flelds (Ohno et al. 2003) and relativistic 
magnetized plasma (Cooray, Melchiorri & Silk 2002). 

The main purpose of this paper is to address the polarization signals discussed above in the context of their statistical 
signature, i.e. the power spectrum, especially in light of the anticipated detection of lensing-induced B-modc signal by Planck 
and other experiments. Our calculated power spectra include the primary Poissonian contribution; for simplicity we ignore 
the smaller contribution due to angular correlation between clusters (Komatsu & Kitayama 1999). The latter peaks on larger 
angular scales that reflect the correlation distance between clusters. 

In section 2 we briefly review the basics of cluster-induced polarization. Calculations of the power spectra are outlined in 
Section 3, followed by results in Section 4, and a brief discussion in Section 5. 



2 POLARIZATION INDUCED BY SCATTERING IN CLUSTERS 

Scattering of the CMB by IC gas changes the radiation temperature along lines of sights to the cluster. In the non-relativisitc 
limit (to lowest order in gas temperature) the thermal SZ effect (Sunyaev & Zeldovich 1972) constitutes a fractional temper- 
ature change 

AT 

— = y9{x) (1) 
g{x) = a;coth(a;/2) — 4 
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(TTrie-^^dl (2) 



where x = hv/{kT) is the dimensionless frequency, Ue and Te are the electron number density and temperature, y is the 
comptonization parameter, ctt is the Thomson cross section, and the integration is along the line of sight. Hereafter we use 
Q for the dimensionless gas temperature kTe/{meC^). 

The second, closely related and smaller kinematic SZ effect is proportional to the line of sight (los) velocity of the cluster, 
Vr, and is independent of frequency, 



AT 

— = - I aTn^Prdl 

Pr = ^. (3) 

C 

Compton scattering can polarize incident radiation if it has a quadrupole moment. The CMB has a global quadrupole 
moment and a non-vanishing quadrupole moment is induced by scattering in the cluster. The degree of linear polarization 
and its orientation are determined by the two Stokes parameters 

Q = ^ J n,dl j sin^ecos2(l>T{e,<t>)dn 

U = ^ j "edZ j sm^esm24>T{e,(l))dQ, (4) 

where 6 and <f) define the relative directions of the incoming and outgoing photons, dQ is an element of integration over the 
solid angle and T{9, 4>) is the temperature of the incident radiation; we use temperature-equivalent units. Since the los is taken 
to be along the z-axis for convenience, the angles 6 and (f> are actually defined with respect to the outgoing photon in this 
system. The average electric field defines the polarization plane with a direction given by 

«=-tan - (5) 
and the total polarization (which is the quantity of interest to us here) is defined as 



P = y/Q2 + U^. (6) 

Another relevant process is the E-B mixing by gravitational lensing. This effect can convert parity-even (E-mode) polar- 
ization to parity-odd (B-mode) polarization (Zaldarriaga 2001). In fact, the largest B-mode signal induced by galaxy clusters 
is due to the gravitational lensing of the primary E-mode polarization. Since CMB lensing by the LSS was extensively studied 
in the past decade, and the power spectrum of CMB polarization due to lensing is readily obtained with Boltzmann codes, 
we do not elaborate on this conversion here (Zaldaxriaga & Seljak 1998). 

In the following subsection we describe the polarization generated by Compton scattering when a quadrupole moment is 
induced by electrons moving either at the cluster peculiar veocity or thermally. When the global quadrupole moment of the 
CMB is taken explicitly into account, then polarization induced by scattering in the gas is treated in the limit when the gas 
is viewed as 'cold', i.e. the second order correction due to random electron motion is ignored. 



2.1 Polarization of EYiUy Isotropic Incident Radiation 

Scattering of the CMB in a cluster at rest in the CMB frame results in local anisotropy due to the different pathlengths of 
photons arriving from various directions to a given point. This anisotropy provides the requisite quadrupole moment; second 
scatterings then polarize the radiation. If the cluster is not resolved, no net polarization would be measured. Nonetheless, it 
is useful to explore the signal associated with double scattering since it is expected to dominate over the other polarization 
signals in rich clusters (Shimon, Rephaeli, O'Shea & Norman 2006) for which is not negligibly small. 

The CMB appears anisotropic in the frame of a non-radially moving cluster; scattering by IC electrons then polarizes it. 
Two polarization components are induced; the first is linear in the cluster velocity component transverse to the los, vt = Btc, 
but quadratic in t; the second is linear in t but quadratic in /3t. The spatial patterns of the various polarization components 
can be readily determined when the gas distribution is spehrically symmetric. The polarization patterns arising from scattering 
off thermal electrons are isotropic in a spherical cluster (Figure 1, left panel) while the corresponding patterns of the kinematic 
components are clearly anisotropic due to the asymmetry introduced by the direction of the cluster motion (Figure 1, right 
panel). More realistic gas distributions, substructure and high internal velocities result in complicated polarization patterns 
(Lavaux, Diego, Mathis & Silk 2004, Shimon, Rephaeli, O'Shea & Norman 2006). 



4 M. Shimon, Y. Rephaeli, S. Sadeh and B. Keating 



T (e=0) Polarization Pattern 



Polarization Pattern 




Figure 1. SZ polarization patterns of the double scattering t^0 (e = 0, x > xq) and t/3j (cluster moving south-west to north-east) 
effects . 



2.1.1 Ther^Pt and T(3f Components 

The degree of polarization induced by double scattering in a non-radially moving cluster was determined by Sunyaev & 
Zeldovich (1980) in the simple case of uniform gas density, 

P=i,r%- (7) 

This polarization component is frequency-independent in equivalent temperature units since it is a first order Doppler shift. 
A more complete calculation of this and the other polarization components was formulated by Sazonov & Sunyaev (1999). 
Viewed along a direction n — {9, (f>), the temperature anisotropy at a point {X, Y, Z), AT{X, Y, Z, 9, 0), leads to polarization 
upon second scattering. The Stokes parameters are calculated from Equation (4), 



Q(X,Y)±iU(X,Y) 



16^ 

3(JT 



dZne{X,Y,Z) I dnsin{9)e^^'''AT{X,Y,Z,9 



= I dZne{X,Y,Z) I dQ.Y2 [9,(i>)AT{X,Y,Z, 



Y^ are the I — 2 and m = ±2 spherical harmonics. The temperature change resulting from first scatterings is 



AT(X, y, Z, 61, < 
T 



OT / dl(X' X ,Z' ,9,4,)n,{X' ,Y\Z\e,<i>)n- P{X\Y' ,Z'), 



and the optical depth through the point (X, F, Z) in the direction (9, 0) is 
r(X,F,Z,6»,0) = aT / n^{X' ,Y' ,Z')dl{X' ,Y' ,Z\9,<I,). 



(8) 



(9) 



(10) 



Q(X, y) and U(X,Y') fully describe the linear 2D polarization field. 

To study the polarization profile we define the following polar coordinates 

X = r cos V' 
Y = r sin -0 

where r is the radial distance from the cluster center. A photon travels a distance 1 before scattering 

X' — X — I sin 9 cos 
y'-y = lsm9sm(i> 
Z'-Z = ZcosS. 

In these coordinates 

Ar(X,y,Z,e,^) ( dl{X',Y',Z',9,<l,)n4X',Y',Z',9,^MX',Y',Z'). 



(11) 



(12) 



(13) 
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Assuming the gas velocity is constant (i.e., internal bulk velocities are relatively low) and that the gas density has the familiar 
/3 profile with index f3 = 2/3, such that the electron density is n(r) = no/{l + r^), where the radial distance from the center 
is now taken to be in units of the core radius, re, we obtain 

AT{x,Y,z,e,ci>) r dl 



^ 7,^0 l + [r2 + «2+2Z(rsin6>cos(V'-<^)+«cos6i)+P] 

where both I, and z are given in units of rc- Writing 
A^^l + r^ + zl 

B = r8in6co8{ip — 4>) + zcosO (15) 
and integrating over I yields 



AT{X,Y,Z,e,4>) « 



TT _1 / B 



(16) 



2 y^/A2—B2 

Inserting this into Eq.(8) and carrying out the 3D numerical integration over z, 6 and 4> we obtain Q and U as functions of r 
and only. The total polarization is then obtained from Eq. (6) P = y^QM-T/^, and by virtue of the fact that P is a scalar 
quantity, independent of ip. By fitting P{r) over a large range of values r (e.g. < ^ < 10) we obtain the following simple 

expression 

PW = ^(n°..M^/3(^^^) (17) 
where the constants a, b and c are obtained by fitting to the numerical 3D integrations described above 
a = 0.57938006 
b = -0.85163944 

c = 0.58122608. (18) 

The second kinematic polarization component is oc r/3f; this component is generated by virtue of the fact that the 
radiation appears anisotropic in the electron frame if the electron motion has a nonvanishing transversal component. From 
the definition of the Stokes parameters Eq.(4) and the fact that the Doppler effect depends on the angle between the photon 
direction and the velocity vector, and by choosing the Z axis to coincide with the direction of the electron velocity, one obtains 
(Chandrasekhar 1950). 

S-^(l-M?)P.(.o), (19) 
where ^2(^0) is the second Legendre polynomial, the expression for the Q parameter of the scattered radiation is 

Q(m) = |t(1 - i4) I P2(MlJ^'o)dii'o ■ (20) 



Here, /io = cos6[), /iq = cos^o, are the cosines of the angles between the electron velocity and the incoming and outgoing 
photons, respectively. Expanding the apparent angular distribution of the radiation in Legendre polynomials, and keeping 
terms up to the quadrupole moment is determined. The polarization of the singly scattered radiation is then calculated 
(Sazonov & Sunyaev 1999) by using Equation (4). The level of this polarization component (Sunyaev & Zeldovich 1980) is 

^= 20(e--l)2 ^^'- (2^) 

In our chosen frame of reference the Stokes parameter U vanishes due to azimuthal symmetry. Therefore, the total polarization 
amplitude, P, is equal to Q and the polarization is orthogonal to l3t. Rclativistic corrections (to the non-rclativistic expression 
of Sunyaev & Zeldovich 1980) were calculated by Challinor, Ford & Lasenby (2000) and Itoh, Nozawa & Kohyama (2000). 
These corrections generally amount to a ~ 10% reduction in the value of Q, and are therefore neglected in our calculations. 



2.1.2 The T^e Component 

Analogous to the component discussed above, double scattering off electrons moving with random thermal velocities can 
induce polarization that is proportional to t^6. The anisotropy introduced by single scattering is the thermal component of 
the SZ effect with temperature change ATt . Its dependence on frequency is contained in the following analytic approximation 
to the exact relativistic calculation (Itoh, Kohyama & Nozawa 1998, Shimon & Rephaeli 2004) 

ATiX,Y,Z,e,4>) ^^^j^ ^^(^^ I „^(,)0(r)*d/. (22) 
i=i ■' 
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In Equation (22) the integration is along the photon trajectory prior to the second scattering, Fi{x) are spectral functions of 
X (Shimon & Rephaeli 2004). The polarization is obtained by inserting this expression in Eq.(8) and repeating the procedure 
described in section (2.1.1). 

The polarization patterns of these two effects are illustrated in Sunyaev & Zeldovich (1980) and Sazonov & Sunyaev 
(1999); as can be seen from Figure 1, the polarization due to double scattering on the thermal plasma is expected to vanish 
when averaged over a beam due to its circular or radial symmetry. When high resoultion, sub-arcminuto experiments will be 
operational the effect will be measured, if sensitivity reaches the few nK level and foreground removal and systematics can be 
controlled at the required level. 

While typical clusters are not perfectly spherical and some residual signal will survive the beam convolution, this signal 
is expected to be very weak; only second order in the cluster ellipticity O(e^). Typically, the correction due to ellipcity is 
comparable to or smaller than few percent. 



2.2 Polarization of Anisotropic Incident Radiation 

The large scale anisotropy of the CMB includes a cosmological quadrupole moment at the level of « 15fiK, as measured by 
the all-sky surveys of COBE and WMAP. Knowing that the probability for generating polarization by scattering in clusters 
is a fraction ~ of the incident quadrupole, and that r ~ 0.01, we expect the resulting polarization signal to be wl5 nK, 
therefore detection would be extremely challenging, even with next generation experiments. 

A similar but smaller component results from scattering of the radiation with a quadrupole moment generated by 
gravitational waves (from tensor perturbations). Its level is likely to be much smaller; the current limit on the tensor-to-scalar 
ratio isT/S < 0.4, deduced from analysis of CMB measurements. The upper (95%) confidence limit from WMAP alone is 0.43 
(Dunkley et al. 2008) and 0.22, if the analysis includes also SN and BAO (Komatsu et al. 2008). We describe the polarization 
induced by the primordial quadrupole (both scalar and tensor contributions) below. The main difference between the two is 
that while density waves do grow during cosmic evolution, tensor perturbations do not, and since we sum the contributions 
from the entire population of clusters, the redshift evolution of the population has to be known reasonably accurately. 



2.2.1 Scalar Contribution 

The global quadrupole moment is imprinted on the CMB by the Sachs- Wolfe (SW) effect at the surface of last scattering, 
and is further boosted by the evolving gravitational fields of density inhomogeneitics - the integrated Sachs- Wolfe (ISW) 
effect. There is no ISW contribution in a purely matter dominated universe, but around recombination the expansion was still 
partially driven by the residual radiation (which constituted about 1% of the energy budget), and more recently (at z < 0.3) 
the expansion became dominated by dark energy, which again caused decay of the gravitational potential. Scattering by IC 
electrons polarizes the radiation at a level proportional to the product of the rms of the primary quadrupole moment and the 
cluster scattering optical depth (Sazonov & Sunyaev 1999). Using the WMAP normalization of the CMB quadrupole moment 
(Bennett et al. 2003), the maximal polarized signal is expected to be ~ 2.6t fiK, and its all-sky average is ~ 60% of this 
value (Sazonov & Sunyaev 1999). 

It has been noted that the dependence of this polarization component on the CMB quadrupole moment could possibly be 
used to reduce cosmic variance (Kamionkowski & Loeb 1997), but this seems doubtful (Portsmouth 2004). Another suggestion 
is that the dependence can probe dark energy models through the redshift evolution of the quadrupole. We note that the 
polarization induced by the CMB quadrupole can be distinguished from other cluster-induced polarization components by 
virtue of its large-scale distribution, reflecting that of the primary CMB quadrupole (e.g Baumann & Cooray 2003), and the 
fact that it is independent of frequency (when expressed in temperature units). The usually quoted value of the quadrupole 
moment Qrms * IS/LtJT , refers to the primordial quadrupole (generated at recombination) but this quadrupole is directly 
afi'ected by the evolving gravitational potential, and since perturbations in this potential can possibly alter the quadrupole, 
it has a redshift dependence imprinted by the evolving potentials. For a power law primordial power spectrum of density 
perturbations P{k) oc fe", the rms value of the quadrupole is (e.g. Hu 2000) 

Qlms{z) = ^5%{1 + zfD\z)Q,i,{dAHoY''Ts^^,{n) 
r M - o /- r[(3-n)/2]r[(3 + n)/2] 

r.w(n) = 3V^r[(4-n)/2]r[(9-n)/2] ^'^^ 
where Sh = 4.2 x 10~®, the growth function is 
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and the Hubble function is 



H{z)/Ho = y^aJT+zj^+oI, (25) 

(where Ho is the present Hubble constant). In the case of a flat power spectrum rsw(«' = 1) = 1 where n is the spectral index 
of scalar metric perturbations, r(a;) is the Gamma function, Qrn and Q,a are the matter and vacuum density today in units of 
critical density, and cIa is the angular diameter distance to redshift z. The polarization level in temperature units is (Sazonov 

& Sunyaov 1999) 

P = ^rQl'^^'r- (26) 



2.2.2 Tensor Contribution 

The scalar metric perturbations dominate the temperature anisotropy at least on small scales but there is also an observational 
upper limit on the tensorial contribution on large scales {I < 100), as well as theoretically expected upper limits (i.e. the latter 
is r < 0.3 for the simplest inflationary models). It is sufficient to gauge this small tensor contribution by the tensor-to-scalar 
ratio r = T/S. The overall normalization should come from observations, and indeed it turns out that the fractional scalar 
perturbations are of order ~ 10^^. Current upper limits from CMB temperature anisotropy experiments should be considered 
weak upper limits, for if inflation indeed occurred following a phase transition in the grand unification (GUT) era, a much 
lower value, r « 0.01, would be expected. It should be noted that, dtie to lensing of the CMB by the LSS,which peaks at 
/ « 1000 but still leaks to the angular degree scales (where the primordial tensor perturbations peak) there is a lower limit to 
the level of tensor-to-scalar ratio (r m 0.001) which can be inferred even from ideal measurements (with noise-free detectors) 
of the primordial B-mode. However, inflation could have taken place before or after the GUT era. The ultimate test will be 
measurements of B-mode polarization, or a direct detection of the stochastic gravitational waves generated during inflation 
(presently considered extremely unlikely). 

As we show below, similar to the polarization induced by the scalar SW and ISW effects, we expect a small polarization 
signal to be induced by scattering, and since the effect scales as Cf oc -^^j^tcnsor^ ^^.^j ^j^^ current upper limit is r < 0.3, 
we expect this signal to be at most ~ 30% of its scalar counterpart. This component contains both E and B modes (due to 
the nature of gravitational waves) and has the potential to assist in tightening the limits on the amplitude of gravitational 
waves. Also, as mentioned above, unlike the scalar quadrupole, the quadrupole induced by gravitational waves does not evolve 
significantly with redshift, i.e. the ISW effect for the tensor modes effectively vanishes. It is sensitive only to the anisotropic 
component of the energy-momentum tensor, and in fact decays by about 10% due to neutrino streaming on quadrupole scales 
(Weinberg 2004), which is a small effect, ignored here. This yields P = ^Ql^i°'^ as in the case of polarization induced by 
Compton scattering of the scalar quadrupole moment. Again, Q^Sf"'^ is obtained from Qrms'"^ (Eq. 23) by setting n — > n — 1 
and the growth function, which depends on redshift, is set to 1, followed by multiplying the resulting power spectrum by r, 
the current upper limit on the tensor-to-scalar ratio. 



3 POWER SPECTRUM 

Although the temperature anisotropy and polarization induced by galaxy clusters are intrinsically non-gaussian and their 
angular power spectra do not contain all the statistical information, these spectra still provide important information on the 
level of the anisotropy and its characteristic scales. The angular power spectrum is essentially the projection of the processed 
3D power spectrum on the 2D celestial sphere. There are two characteristic angular scales in the problem; the angular size of a 
cluster and the typical angular correlation angle on the sky between neighboring (sufficiently rich) clusters. The first is typically 
a few arcminutcs aud the correlation angle is ~ 1° (if only rich clusters are considered) which correspond to I « 2000 — 3000 
and / « 200, respectively. Their relative importance depends on the corresponding 3D power spectra which consists of the 
Poisson and correlation (clustering) terms. The contribution due to correlations is typically an order of magnitude smaller 
(e.g., Cooray, Baumann & Sigurdson 2005) and is not considered here. The power spectrum of the anisotropy due to the 
population of clusters is 

where ^ is the Fourier transform of the polarization generated in each cluster 

5(1) = j d^e^ie)e''-^. (28) 
The above integral is over the mass function of clusters, for which the Press-Schechter distribution (or one of its variants) 
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is usually adopted 
where 

F(M) = ./fe-^-f, (30) 
M{R) = ^pbR^ and the mass variance is 

a'{M;z) = D'{z) J !^t!£^\w{kR)\\ (31) 

Here n = Sc{z)/a{z) is the critical overdensity for collapse in terms of the mgiss variance at redshift z, and pt, is the background 
density at 2; = 0, and W{q) = -^[sm{q) — qcos{q)] is a top-hat window function. The processed density fluctuation power 
spectrum P{k) = Ak"T^{k) is given in terms of the transfer function T(k) (which is specified below). The normalization 
constant is determined in terms of the observationally deduced value of as, the mass variance on a scale of 8/i~^ Mpc. 
The mass variance evolution with redshift is given by 

-(^;^) = ^B4i4^ (32) 

^ ' g[Q.mm 1 + z ^ ' 

where 

^^{z) ^ n^(Q){l + zf . (33) 
For the function (;[f2,„(2;)] we use the approximate expression (Carroll, Press & Turner 1992) 

a\0 (zW = 2.5n,n(g) ,„ 

>^ n^(0)4/7_nA + (i + Om(;?)/2)(i + nA/7o)' ^ ' 

We adopt the standard CDM transfer function 

nk) = ^"^^+3^^^^^^ [1 + 3.899 + (16.1g)^ + (5.46g)^ + {Q.llqf]-''' (35) 

with q = k/{Q.mh^)Mpc~'^ (Bardeen et al. 1986). 

The virial relation is used for the gas temperature is obtained assuming virialization 

3f3r,{M;z) 

where the core-radius is obtained from the spherical collapse model 

1/3 

r-ciM; z) ° 



1 + z 



(37) 



with ro = 1.69 Mpc/(/ip) where p is the ratio between the virial to core radius Rv/rc (taken here to be 10). The distribution 
of cluster bulk velocities is derived from the continuity equation 

vL. = J ^Pik). (38) 

The calculation of ^(1) is significantly simplified by focusing on the total polarization which unlike Q and [/ - is a scalar 
field (i.e. a function of r only, as in Eq. 17 or the 2D-projected /3-profile), and - as in the case of temperature anisotropy - we 
can employ the expansion of scalar-valued 2D plane wave in cylindrical Bessel functions 

e*'-* = ^i'j((/6l)e'('^-'^'' (39) 

I 

to obtain 

1(1) = 27ry" edejo{w)i{e) (40) 

where ^(6') is the total polarization. Also, use of a /3-profile for the density implies ^{9) — ^{9/6c), with 6c having explicit M 
and z dependence 6c = 6c{M; z) through the spherical collapse model. Now the mass and the radius are related through the 
density which is fixed to be l^iit^ ph{z), where ph is the background density following the simple spherical collapse model. This 
density depends on the redshift and therefore the electron density depends on r, M and z. 

To obtain ^(6) we simply take the square root of Q and U added in quadrature at an angular distance 9 = r/dA from the 
cluster, where r and dA are the physical distance from the cluster center (r = \/X^"+Y^), and the angular diameter distance 
to the cluster, respectively. 



Power Spectra of CMB Polarization by Scattering in Clusters 9 



Cluster Induced Polarization (143 GHz} 



Cluster Induced Polarization (217 GHz] 




Cluster Induced Polarization (353 GHz) 



Cluster Induced Polarization (545 GHz) 



- Prim. Quad. S 

- Prim. Quad. T 



\ 




- Prim. Quad. S 

- Prim. Quad. T 



\ 

v. \ 




Figure 2. SZ polarization at Planck HFI frequencies; Dot-dashed red lines are the B-modes from the primordial tensor perturbations 
with tensor-to-scalar-ratios r = 0.1, 0.01 and 0.001 respectively. Dot-dashed blue curve is the E-mode polarization. Dot-dashed green 
curve is the B-mode from lensing conversion of the E-mode by the LSS. Assuming that this signal can be cleaned to 1 part in 40 we plot 
the residual B-mode from lensing (dashed green curve). Also shown are five SZ (total-) polarization signals: primordial scalar quadrupole 
(blue), upper limit (r=0.3) for the effect due to primordial tensor quadrupole (black), r^© (green), (red) and (magenta). 



4 RESULTS 

All power spectra presented in this work were calculated assuming the ACDM with the WMAP5 best-fit parameters (Komatsu 
et al. 2008): baryon, dark matter and dark energy densities (in units of the critical density) Q^h^ = 0.0227, Qch^ ~ 0.1099, 
f^A ~ 0.742, respectively, and scalar spectral index, Hubble constant in units of 100 km/sec/Mpc — 0.963, h = 0.72, 
respectively. The mass variance parameter is ag — 0.8. 

Power spectra of the cluster polarization components are shown in Fig. 2 together with the corresponding ones of the 
primary anisotropy. The largest CMB polarized power is the primary E-mode polarization shown by the dot-dashed blue line 
in Fig. 2. Lensing of this E-mode anisotropy by the LSS produces B-mode pattern (dot-dashed green line in Figure 2) on a 
level which severely contaminates the primordial B-mode signal from inflation, shown by the dot-dash red lines in Figure 2 for 
three values of the tensor-to-scalar ratio r. As of yet the only information we have on r, which determines the amplitude of the 
B-mode from primordial gravitational waves, is an upper bound of r < 0.3. Various techniques have been devised to remove 
the contamination due to the LSS from the primordial signal to allow the determination of the energy scale of inflation (e.g. 
Hu & Okamoto 2002). However, below r ~ 0.001, the residual lensing still dominates the B-mode power even after the lensing 
contribution is optimally reduced by a projected factor of ~ 40, as demonstarated by the dashed green curve in Figure 2. 
Furthermore, for r ^ 0.01 the detection of the primordial B-mode is unlikely due to foreground contamination. Accordingly, 
the gravitational wave signals for r = 0.1, r = 0.01 and r — 0.001 are shown for comparison. 
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Our results for the total polarization power spectra induced by scattering in clusters are shown by the solid lines in 
Fig. 2. Both the polarization due to scattering of the primordial quadrupole (scalar and tensor), the t^/Jt effect, and the 
E/B conversion arc independent of frequency, and thus are spectrally indistinguishable from the primordial signal. Separation 
between these contributions can possibly be based on their different statistical properties (e.g. the characteristic polarization 
pattern of the former, and non-gaussianity in the later case). Clearly, the largest B-mode signal comes from the E/B mixing 
by lonsing. The t0^ and the t^Q components do depend on frequency (as discussed in the previous section) and can be 
removed from the primary CMB signals by multifrequency observations. 

All our results are based on the ACDM with WMAP5 parameters. However, since the SZ effect is a sensitive function of 
matter clustering, as, and since high-resolution CMB experiments, such as BIMA, CBI and ACBAR, which are sensitive to 
cluster angular scales infer higher us than the WMAP5 value of ~ 0.82. With as = 0.9, as suggested by galaxy surveys, the 
cluster polarization power levels are higher by a factor of at least ~ 2 (see e.g. Sadeh, Rephaeli & Silk 2007). Polarization 
from the tensorial quadrupole (black solid line in Fig. 2) is only an upper bound corresponding to r = 0.3 (close to the current 
upper limit 0.22). 

5 DISCUSSION 

All cluster components calculated here are sub- to few nK and are smaller than the B-mode signal from lensing by 3 — 4 
orders of magnitude, too weak to have any noticeable impact on parameter estimation, even those which depend on the 
largo scale structure, e.g. neutrino mass Mi,, dark energy equation of state, w, etc. However, assuming that the lensing signal 
can be removed to 1 part in 40 by invoking optimal filters based on higher order statistics, cluster signals could potentially 
contaminate the residual lensing-induced B-modes signal on the few-percent level, again - too weak to have an impact on 
inflationary models. 

Polarization signals from individual rich clusters may be detectable. While these signals may be negligible when averaged 
out on the full sky, current ground-based endeavors to detect the B-mode polarization typically target only few percent of the 
sky. Although these 'radio-quiet' patches arc optimally chosen to minimize the noise, it is not entirely excluded that CMB 
polarization by individual clusters will have an overall effect of more than the conservative few percent residual contamination 
shown in Figure 2. Also, it is important to reiterate here that the method of Hu & Okamoto, mentioned above in the context 
of dolcnsing the sky from the LSS-inducod B-mode, indeed rests on the assumption that the unlensod signal is gaussian and 
that the expected small level of non-gaussianity is due to lensing by the LSS. Therefore, it cannot be employed to remove 
the non-gaussian SZ polarization from clusters (which indeed was not calculated in this work). This will require, in principle, 
a customization of similar spatial-filtering methods to the case of galaxy clusters; a problem we have not addressed in this 
work, mainly because it is expected to be negligibly small for polarization. 

The morphology of IC gas typically shows a finite level of ellipticity. It is of interest to assess the implications of a small 
level of ellipticity on the small-scale power spectrum of the r^G polarization component. As Figure 1 illustrates, the effect of 
double scattering by IC gas in a spherically symmetric cluster will result in vanishing polarization at the cluster center, due to 
the fact that no quadrupole is generated there from the first scattering. As a result, the r^O component drops on very large 
I. However, cluster ellipticity changes this behavior because first scattering of photons travelling along the major and minor 
axes will result in a few percent quadrupolar change in their corresponding SZ-temperatures, and upon second scattering this 
small quadrupole will be polarized even at the cluster center. The leading order correction to the r^O power spectrum will 
be O(e^), rather than 0(e) (since both positive and negative ellipticity statistically avcrgac out to (e) = 0). In addition, the 
main qualitative effect will be to boost power at the cluster center as explained above, but this will probably be beyond the 
reach of even the highest resolution next generation experiments. 

ACKNOWLEDMENTS 

We acknowledge using CAMB to calculate the primordial power spectra. BK gratefully acknowledges support from NSF 
PECASE Award AST-0548262. 

REFERENCES 

Audit, E., & Simmons, J. F. L. 1999, MNRAS, 305, L27 

Baumann, D., & Cooray, A. 2003, Now Astronomy Review, 47, 839 
Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 
Bennett, C, et al. 2003, ApJS, 148, 1 

Carroll, S. M., Press, W. H., & Turner, E. L. 1992, ARAA, 30, 499 



Power Spectra of CMB Polarization by Scattering in Clusters 



Challinor, A. D., Ford, M. T., & Lasenby, A. N. 2000, MNRAS, 312, 159 
Chandrasekhar, S. 1950, Oxford, Clarendon Press, 1950 
Cooray, A., Melchiorri, A., Silk J. astro - ph/0205214 

Cooray, A., Baumann, D., & Sigurdson, K. 2005, Background Microwave Radiation and Intracluster Cosmology, 309 

Dunkley, J., et al. 2008. [arXr v:0803.0586 

Gibilisco, M. 1997, Ap & SS, 249, 189 

Hinshaw, C, et al. 2008, ,arXiv:080_3_-0I32 , 

Hu, W. 2000, ApJ, 529, 12 

Hu, W., & Okamoto, T. 2002, ApJ, 574, 566 

Hu, W., DeDeo, S., & Vale, C. 2007, ArXiv Astrophysics e-prints, [MXiv:astro-ph/0701276| 

Itoh, N., Kohyama, Y., & Nozawa, S. 1998, ApJ, 502, 7 

Itoh, N., Nozawa, S., & Kohyama, Y. 2000, ApJ, 533, 588 

Kamionkowski, M. & Loeb, A. 1997, PRD, 56, 4511 

Komatsu, E., & Kitayama, T. 1999, ApJL, 526, LI 

Komatsu, E., et al. 2008, arXiv:0803.0547 

Lavaux, C, Diego, J. M., Mathis, H., & Silk, J. 2004, MNRAS, 347, 729 
Lewis, A., & King, L. 2006, PRD, 73, 063006 

Ohno, H., Takada, M., Dolag, K., Bartelmann, M., & Sugiyama, N. 2003, ApJ, 584, 599 

Portsmouth, J. 2004, PRD, 70, 063504 

Sadeh, S., Rephaeh, Y., & Silk, J. 2007, MNRAS, 380, 637 

Sazonov, S. Y. & Sunyaev, R. A. 1999, MNRAS, 310, 765 

Seljak, U., & Zaldarriaga, M. 2000, ApJ, 538, 57 

Shimon, M., & Rephaeh, Y. 2004, New Astronomy, 9, 69 

Shimon, M., Rephaeh, Y., O'Shea, B. W., & Norman, M. L. 2006, MNRAS, 368, 511 
Silk, J. 1968, ApJ, 151, 459 

Sunyaev, R. A. & Zeldovich, 1. B. 1972, Comm. Ap. Sp. Phys., 4, 173 

Sunyaev, R. A. & Zeldovich, 1. B. 1980, MNRAS, 190, 413 

Weinberg, S. 2004, PRD, 69, 023503 

Zaldarriaga, M., & Seljak, U. 1998, PRD, 58, 023003 

Zaldarriaga, M. 2001, PRD, 64, 103001 



